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In this paper the theoretical model is built for ZEpHyR (ZARM Experimental Hybrid Rocket) main engine which 


is being developed at ZARM institute, Bremen, Germany. The theoretical model is used to estimate the tempera- 


ture of exhaust gas. The Conjugate Gradient Method (CGM) with Adjoint Problem for Function Estimation itera- 


tive technique is used to solve the Inverse Heat Conduction Problem (IHCP) to estimate the heat flux and internal 


wall temperature at the throat section of the nozzle. Bartz equation is used to calculate the convective heat trans- 


fer coefficient. The exhaust gas temperature is determined using the estimated heat flux, the wall temperature at 


internal surface of nozzle and the heat transfer coefficient. The accuracy of CGM iterative scheme to solve the 


IHCP is also investigated and its results are presented. 
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Introduction 


In the past couple of decades, development of the hy- 
brid rockets has increased. The major reason being hybr- 
id rockets are much safer and easier to use than liquid 
and solid rockets in terms of complexity and controllabil- 
ity. At ZARM Institute — Univeristét Bremen, a hybrid 
rocket engine is being developed under the framework of 
the DLR-STERN (Studentische Experimentalrakten) 
program by DLR, German Aerospace Center. The ZE- 
pHyR (ZARM Experimental Hybrid Rocket) engine runs 
on paraffin as solid fuel and liquid oxygen as the oxidiser 
to power it to an altitude exceeding 12 km [6]. One of the 
factors significant to the success of the rocket engine 
depends on the nozzle. The nozzle has to survive extreme 
temperature and pressure, where the combustion products 
(exhaust gas temperature) entering the nozzle can reach 
up to 3500K and a pressure of up to 35 bar. Only handful 
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of materials are present that can withstand these condi- 
tions without active cooling. Knowledge of the total 
heat transferred by the combustion gases is essential in 
order to find the temperature distribution in a rocket noz- 
zle. Estimation of heat flux on nozzle wall is important to 
achieve best possible thermal protection system. In this 
paper, a hybrid approach is developed to estimate the 
temperature of the exhaust gas of ZEpHyR engine. Here, 
the hybrid approach is combination of both an iterative 
scheme and Bartz equation to estimate the exhaust gas 
temperature. 


Problem Description 


It is very difficult to determine the unknown boundary 
condition by direct measurements in a rocket nozzle as it 
is subjected to very high temperatures. Due to extreme 
environmental conditions many experimental problems 


Nomenclature 


T; Inner radius [m] Pr Prandtl number 

f Outer radius [m] M Mach number 

Ar Node thickness [im] y Specific heat ratio 

D, Throat Diameter [m] J Number of thermocouples 
Tt Throat radius of curvature [m] L Number of measurements 
At Throat area [m] € Tolerance 

Ac Chamber area [m] (o Standard deviation 

At Time difference [s] 0) Random variable 

T, Initial temperature [K] Subscripts 

Y Measured temperature [K] 0,1.. Node identifier 

P, Chamber pressure [bar] Ww Wall surface 

c* Characteristic velocity [ms"'] g Exhaust gas 

p Density [kgm] s Stagnation 

k Thermal conductivity [Wm'K"] f Final value 

h Convective heat transfer coefficient [Wm°K"] e Exact 

q Heat flux [Wm] j Thermocouple number 
Cp Specific heat at constant pressure [kg "K "] 

g Gravitational acceleration [ms] 

u Viscosity [Pas] 


arise in implanting a sensor at the internal wall surface of 


the nozzle, to take heat flux or temperature measurements. 


Moreover implanting the probe at the internal surface 
disturbs the boundary condition and the accuracy of the 
measurements. But it is easier to accurately measure the 
temperature at interior location or outer surface of the 
rocket nozzle. This type of problem is called inverse 
problem or Inverse Heat Conduction Problem (IHCP) 
and requires an iterative technique for its solution. By 
using inverse technique, unknown boundary condition 
can be estimated indirectly from the measured tempera- 
ture within the nozzle. In this paper one dimensional 
transient heat transfer mathematical model which de- 
scribes the processes involved is established to solve the 
IHCP in order to estimate the internal wall temperature 
and heat flux. 


The mathematical model 


In general, the mathematical theoretical model repre- 
sents the physical system and processes involved in it 
and also it should be able to accurately predict the system 
behaviour. In this research the one dimensional theoreti- 
cal model of the nozzle system is built using the heat 
transfer theory. Transient heat transfer model is built by 
taking the nozzle throat as hollow cylinder of finite thic- 
kness where it is convectively heated at the internal sur- 
face and perfectly insulated at the outer surface. 

Consider a hollow cylinder of finite thickness with in- 
ternal radius ‘r; and external radius ‘r,’ and assume con- 


stant thermal properties. Cylinder is perfectly insulated at 
T= To. A heat flux ‘q (t)’ is subjected at r = r;, at a time t > 
0. Figure 1 shows the geometry of the model built. 
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Fig. 1 The model geometry 


The mathematical formulation for the one dimensional 
problem considered is given in Equations (1.1) to (1.4). 
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The implicit method (backward difference) is used to 
numerically represent the above system in finite differ- 
ence form. Implicit method uses both the current and 


future state of the system to solve the equation. The tran- 
sient heat transfer equations for the model are written in 
finite difference form for ‘n’? number of nodes, current 
time ‘i’ and future time ‘i+1’. The equations are repre- 
sented in tri-diagonal form as follows. 
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Bartz equation: 

Prior knowledge of the heat transfer coefficient is re- 
quired to estimate the exhaust gas temperature. To deter- 
mine the heat transfer coefficient Bartz equation [[1]] is 
used and is as follows. 
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The values of Prandtl number, viscosity, specific heat 
and characteristic velocity are obtained from NASA 
Chemical Equilibrium with Applications (CEA) program 
[5]. 

Using the relationship given in Equation (4) the ex- 
haust gas temperature, T, can be calculated. 

q(r,t) =h(T, -T,,) (4) 


The ‘q(r,0)’ and ‘Tw’ values in above equation are 


calculated from IHCP while the heat transfer coefficient 
‘h’ is calculated from the Bartz equation. 

Using multiple thermocouples the transient tempera- 
tures can be measured at outer surface and at various 
depths of the nozzle wall. The IHCP is solved using the 
iterative scheme called "Conjugate Gradient Method with 
adjoint problem for function estimation" [2] to estimate 
the internal wall surface temperature and heat flux, using 
the temperature history measured at interior locations and 
outer surface of the rocket nozzle. This iterative tech- 
nique is used to solve for the heat flux as no prior infor- 


mation about its functional form is available. This me- 
thod is applied to the mathematical model developed. 
The iterative scheme used here comes under the class of 
Alfanov’s iterative regularization methods [3].The con- 
vective heat transfer coefficient is calculated using the 
Bartz equation. Using the results of the iterative method 
and Bartz equation, the exhaust gas temperature is calcu- 
lated. 


Results and Discussion 


The rocket nozzle of ZEpHyR is made of Molybde- 
num having r; = 0.0103m and r, = 0.0305m. The proper- 
ties of Molybdenum are summarized in Table 1. Five 
thermocouples are used to take the temperature mea- 
surements, one at outer surface and four thermocouples 
at varying depths of the nozzle wall as shown in Figure 2. 
The parameters of ZEpHyR are used to estimate the heat 
flux and nozzle internal wall surface temperature using 
CGM for n = 21 nodes. The results obtained along with 
the data obtained from NASA CEA Program are used to 
solve Bartz equation to approximate the heat transfer 
coefficient. Using all these values the exhaust gas tem- 
perature is computed. 


Table 1 Properties of Molybdenum [4] 
Property Value 
k (Wm'K’") 138 
cp kg! K») 250 
p (kgm) 10220 


all dimensions in mm 


Fig. 2 Thermocouples placed at varying depths at nozzle 
throat 


Verification for accuracy of CGM method: 

To demonstrate the accuracy of CGM method in esti- 
mating the heat flux and surface wall temperature at 
boundary, two test heat flux functions are assumed to 
collect the exact temperatures’ (Te) data by the direct 
method given by Equations (1.1) to (1.4) for the total 
duration of 60s. 

The thermocouples used to obtain temperature mea- 


surements during rocket nozzle testing have measure- 
ment error. But the exact temperatures obtained by using 
the test fluxes do not include these errors. Therefore, 
random measurement errors are introduced in the exact 
temperatures obtained by direct method. This is done by 
Equation (5) [2]. 
Y=T,+ @o (5) 
In Equation (5), ‘o° is the standard deviation and ‘œ’ is 
the random variable having normal distribution with zero 
mean and unitary standard deviation. The random varia- 
ble ‘œ’ is in the range -2.576 < œ < +2.576 for the confi- 
dence bound of 98.9%. Y is the measured temperature 
with errors [2]. The measurement errors considered are 
‘o’ of 3 K, 5 K and increased to 10 K. The stopping crite- 
rion is calculated for each standard deviation from Equa- 
tion (6). 
E= Ot (6) 
The relative root means square error for heat flux 
‘ehrms and internal wall temperature ‘e,;’ is calculated to 
find the difference between the estimated results obtained 
from CGM and the actual results from direct method 
using test heat flux functions. The relative means square 
error is obtained by using the Equations (7). 
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In Equations (7), ‘L’ denotes the total number of mea- 
surements, ‘ġ'(r,, 7)’ and ‘7'(7,,/)’ represents heat flux 


and internal wall temperature respectively. 

Test case I (Step heat flux): The step change in the 
heat flux defined in Equation (8) is used to obtain the 
exact temperatures from the direct problem. Random 
measurement errors are introduced having ‘o’ of 3 K, 5K 
and10 K in the exact temperatures using Equation (8). 


0 0<ż<10 
g =45.5x10°10<1<50 (8) 
0 50<1t< 60 


The comparison of results obtained for ‘o’ of 3 K, 5 K 
and10 K is shown in Figure 3, Figure 4 and Figure 5 re- 
spectively. It can be seen that with increase in the stan- 
dard deviation, increase in the difference in results be- 
tween the test and estimated heat flux and wall tempera- 
ture also increases. The root means square errors are cal- 


culated for the three standard deviations and are summa- 
rized in Table 2. 

The general trend observed from the values of Table 2 
is that with increasing ‘o’, the error increases but maxi- 
mum error observed is only 0.0113 %. Study of the re- 
sults shows that the solution of the inverse heat conduc- 
tion problem using CGM iterative scheme remains stable 
with increase in errors. 
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Fig. 3 Comparison of test and estimated step heat flux and 
temperature of internal wall for ‘o’=3 


Test case II (Ramp heat flux): The linear variation in 
test heat flux is used to obtain the exact temperatures 
from direct problem. Random measurement errors are 
introduced in the exact temperatures using Equation (9). 

_ {1x10°¢-1000 0<1<30 
q= : (9) 
-1000002+6x10° 30<1<60 

The comparison of results obtained for ‘o° of 3 K, 5 K 
and 10 K is shown in Figure 6, Figure 7 and Figure 8 
respectively. Analysis of results show that as the standard 
deviation is increased the difference in results between 
the test and estimated heat flux and internal wall temper- 
ature also increases. The root mean square errors are 
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calculated for the three standard deviations and are sum- 
marized in Table 3. 

The general trend observed from the values of Table 3 
is that with increasing ‘o’, the error increases but maxi- 
mum error observed is only 0.0505 %. It can be observed 
that the results of the inverse problem using CGM me- 
thod remain stable as errors are increased. 
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Fig. 4 Comparison of test and estimated step heat flux and 
temperature of internal wall for ‘o’=5 


Table 2 Root mean square error in heat flux and internal wall 
temperature for step heat flux 


Standard deviation ’6’ Ehrms Ctrms 


6=3 0.0011 0.0011 
6=5 0.0018 0.0019 
6=10 0.0113 0.0084 


Estimation of exhaust gas temperature: 
The exhaust gas temperature is estimated using Equa- 


tion (4). The ‘ q(r,t) > and ‘Tw’ in this equation are ob- 
tained from the results of test case II (ramp heat flux) for 
‘o’ of 3. Further, the Bartz equation is used to calculate 


ChinaXive (ERAT 


the heat transfer coefficient ‘h’. 
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Fig. 5 Comparison of test and estimated step heat flux and 
temperature of internal wall for ‘o’=10 
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Fig. 6 Comparison of test and estimated ramp heat flux and 


temperature of internal wall for ‘o’=3 
Table 3 Root mean square error in heat flux and internal wall 
temperature for ramp heat flux 


Standard deviation ‘ó’ Ehrms Ctrms 
6=3 0.0212 0.0062 
6=5 0.0306 0.0097 
6=10 0.0505 0.0182 
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Fig. 7 Comparison of test and estimated ramp heat flux and 
temperature of internal wall for ‘o’=5 


Calculation of heat transfer coefficient: The heat 
transfer coefficient is estimated at the throat section of 
the nozzle by using Equation (3). The parameters re- 
quired for its calculations are obtained from NASA’s 
CEA program for the oxidizer to fuel mixture ratio of 2.1 
and the chamber pressure of 3.5x10° Pa. The oxidizer 
and fuel used in ZEpHyR are liquid oxygen and paraffin 
respectively. The properties and transport parameters 
obtained from NASA’s CEA program are listed in Table 4. 

The heat transfer coefficient is calculated using the 
information given in Table 4 combined with the internal 
wall temperature obtained from CGM iterative scheme. 
In equation (4), heat transfer coefficient is inversely pro- 
portional to the wall temperature and similar trends are 
observed in the plots of heat transfer coefficient shown in 


Fig. 8. 
Exhaust gas temperature: Equation (4) is re-arranged 
as follows to find the exhaust gas temperature. 


q(r,t) 
T= Ea +T, (10) 
The exhaust gas temperature profile over the time of 
60s is shown in Fig. 10. It can be observed that the ex- 
haust gas temperature shows almost linear variation with 
time which is similar to the variation of test ramp heat 


flux applied as shown in Fig. 6. 
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Fig. 8 Comparison of test and estimated ramp heat flux and 
temperature of internal wall for ‘o’=10 


Table 4 Properties and transport parameters of combustion 
gas [6] 


Property Value 
ì (Pas) 9.8873x10°5 
cp (kg! K>» 2028 
D (m) 0.0206 
r (m) 0.0103 
c* (ms') 1768.1 
Pr 0.6196 
Ac/At 1 


a 1.1572 
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Fig. 9 Variation in heat transfer coefficient with time and 
internal wall temperature 
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Fig. 10 Estimation of exhaust gas temperature for ramp heat flux 


Conclusion 


The theoretical model developed for ZEpHyR nozzle 
is used to obtain the temperature of exhaust gas by ap- 
plying Conjugate Gradient Method iterative technique to 
obtain the solution for the inverse heat conduction prob- 
lem and Bartz equation to compute the heat transfer coef- 
ficient respectively. Conjugate gradient numerical itera- 
tive technique used to estimate the heat flux has the ad- 
vantage that it requires no prior knowledge about the 
functional form of boundary heat flux and has high rate 
of convergence. Also high accuracy of results is obtained 
from CGM that shows that it is less sensitive to mea- 
surement errors. The hybrid approach adopted in this 
research can be accurately used to estimate the exhaust 
gas temperature of ZEpHyR. The ZEpHyR nozzle has 
not reached the testing phase yet due to which measured 
data for temperature is not available from the thermo- 
couples. The test heat flux used shows promising results 
and the technique presented here can be applied to esti- 
mate the exhaust gas temperature using the actual test 
thermocouple data. 
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